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NEUTRAL BEAM DOSE AND SPUTTERING CHARACTERISTICS IN 


AN ION IMPLANTATION SYSTEM 

By 

A.S. Roberts, Jr. 1 , R.L. Ash 2 and 
M.H. Berger 3 

1. INTRODUCTION 

Accurate determination of ion beam dose is important in the production 
of ion implanted semiconductor devices. A small but significant percentage 
of neutral atoms in the source beam can lead to inaccuracies in impurity doping 
and localized sputtering of the semiconductor target in a high-vacuum ion 
implantation system. The work reported here discloses a technique and instru- 
ment design for calorimetric detection of the neutral atom content of a 60 
keV argon ion beam. A beam sampling method is used to measure local heat flux 
to a small platinum wire at steady state; integration of power density profiles 
leads to a determination of equivalent neutral beam current. The fast neutral 
production occurs as a result of charge transfer processes in the region of 
the beam system between analyzing magnet and beam stop where the pressure 
remains less than 10" 5 torr. A description of the neutral beam detector is 
given in section 2 along with a presentation of results. 

As a part of this work, performed for the Instruments Techniques Section, 
Measurements Physics Branch of the Instrument Research Division, an elementary 
analysis of sputter material transport from target to substrate was performed; 
the analysis relates to semiconductor sputtering. This work with a cursory 
review of semiconductor sputtering literature is described in section 3. All 
experiments were performed using a high vacuum ion implantation system in the 
Instruments Techniques Section laboratories at NASA Langley Research Center. 
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2. DESIGN OF NEUTRAL BEAM, DETECTOR 

A conventional plasma source was used to generate ionized argon; the 
cathode was tungsten. Extracted tons were formed into a beam via an electro- 
static lens system and the beam was magnetically analyzed to produce a singly 
ionized argon beam in the target region. A target chamber schematic is shown 
in Figure 1. The target chamber is offset upward so that when the ion beam 
is electrostatically deflected upward and focused onto the center of the beam 
stop the neutral beam will continue off-center to the end of the target chamber. 
As indicated, the beam detector is moved up through the neutral beam. Figure 2 
is a photograph of the target region assembly. The target chamber is left of 
center with the beam detector inserting from the bottom; below the large 
central tee a noble gas sputter ion pump is mounted; the beam enters from the 
right through the metal bellows; the vertical offset is visible by observing 
the bellows position. Finally, beam detection instrumentation is seen in the 
foreground. 

2.1 Description of the Wire Detector 

The problem was to devise an instrument to detect an equivalent neutral 
beam current down to -v .1 pa. The neutral beam is formed between the analyzing 
magnet and the vertical deflection plates (a distance of ^ 75 cm); beam energy 
is assumed uniform at the value set for the ion beam, 60 keV.' 

A number of techniques have been employed to detect neutral s,< commonly 
involving deposition of energy or momentum impulse (Huddlestone and Leonard, 
chap. 12) 1 . On the other hand, Berkner, et al. used a pyroelectric ceramic 
detector to determine neutral beam fluxes; however, the method yielded no beam - 
structure and required simultaneously a measurement of ion beam current. A . _ ; 
note by Ormrod disclosed an energy deposition technique where an ion beam was 
scanned with a tungsten wire and the resulting black-body radiation was 
measured. Although the method gave beam structure information it is only use- 
ful for high beam currents (significant heating of wire). 

An energy deposition technique was ultimately used. A 1 mil diameter 
(.0254 mm) platinum wire, ^ 2 cm in length, was swept through the beam, sub- 
tending the beam in a chord-wise fashion. Since beam diameters were 7 to 10rrmm> 

^See Reference list. \ j 




Figure 1. Target chamber schematic. 
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Figure 2. Target region assembly. 





a number of sampling points could be gotten over the beam cross-section. A 
photograph of the wire detector and its mount is shown as Figure 3; the detector 
assembly inserts into the bottom of the target chamber as sketched in Figure 1. 
The neutral beam energy flux can be inferred from the measured average change 
in platinum wire resistance; the neutral beam dose then follows if it is as- 
sumed that (1) neutral atom energy is equivalent to that of the ion beam and 
(2) kinetic energy of neutral atoms is completely absorbed when intercepted by 
the wire. 

2.2 Analysis for Data Reduction 

The platinum wire is heated by the beam and, when steady-state is reached, 
the absolute resistance is measured with an accurate Kelvin Bridge. If the rise 
in average wire temperature, with respect to its measured end temperature, is 
limited to a few hundred degrees centigrade, a linear relation holds between 
average incremental resistance and average incremental wire temperature; 




<r> 

a r ref 


a = .0039 °C _1 (1) 

(r f is resistance of wire 
prior to beam heating) 

where platinum properties are taken from a monograph by Caldwell. Note that 
an accuracy for <T p > determination of ^ 2.5% is implied. Figure 4 outlines 
the experimental technique for determining <T p > from which beam power density 
can be inferred for a given position of the wire in the beam. 


A steady-state heat balance is performed, and an average heat flux can 
be calculated for a given <T p >. The one-dimensional governing equation is 
fourth degree, subject to fixed wire end temperature. 

TT R2 k ^ - 2^Rea(T 4 - TA 4 ) + 2RH = 0 (2) 

Terms have dimensions watts per unit length of wire. Energy loss from the wire 
is attributed to axial conduction (x - direction) and black-body radiation; 
k is thermal conductivity, e emissivity, a Stephan-Bol tzman constant, TA 
enclosure temperature, and other symbols are made clear in Figure 4. Boundary 
conditions are, 


£ <°> - 0 
T(L) = TEND 


(3) 
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Vive deteotov assembly . 
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Figure 4. Method of determining beam power density. 

(a) Platinum wire subtends cord in beam of diameter 
BDIA » 2R/(b) energy of beam, Hj (watts/cm 2 ) is 
absorbed. until losses balance gains, (c) a steady 
state temperature profile obtains with <T p > the 1 
experimental ly.'determined average incremental 
temperature,. Xx , the half-heated length; and -L., ! \~ 

the half-wire length. - .. ! 



The problem of data reduction amounts to the determination of the power 
density of an incident beam given: 

(1) the average temperature of the wire 

(2) the position of the wire in the beam 

(3) the beam diameter 

/ \ . . 

/ (4) the rod end temperature, TEND 

(5) the ambient (enclosure) temperature, TA 
A great deal of information is revealed by studying the solution tc? the linear- 
ized equation. 

- ’^Define 

T(x) = TA (T p (x )/tJA !+ ilj) (4) 

where T p (x)^TA is small compared with unity. By performing the Taylor Series 
expansion of (T(x) 4 - TA 4 ), replacing T(x) with the previous definition and 
applying the perturbation condition a new differential equation in T p (x) results. 




where. 


y" 


H(x) = H ls constant 0 £ x £ Xx 
0, x. 1 < x < L 
The new boundary conditions are. 


(5) 


0 < x < Xx 


dT^O) = 0 
dx 

dT,(*i‘) * £,(xi + ) 

dx dx 


( 6 ) 


y*r> - T (xi + ) 

Xi < x < L | y H 

T (L) = TEND - TA 


"S 


: •' : "an pa ah 

Solutions can be obtained as follows: 

For 0 £ x £ xx 

T p (x ) = (Qx/Bx) + C 2 cosh x) 


(7) 


(8) 
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For 


Xi < x _< L 


Tp(x) = F l sinh (L-x)) + "(TEND - TA) 


( 9 ) 


cosh {/B! (L - x)) 


The incremental temperature distribution is algebraic in Qi (“Hi ) . The average 
incremental temperature for a particular Hi value is obtained by integrating 
equations (8) and (9); 


L 

'V = r j T P (x) dx 

0 

xi L 

= [ | T p(x) dx + f T p (x) dx 

0 xj 


(10) 


<T > 

P 


f 1 (Qi/Bi) x x + — sinhtiWxi) ; 
L I /B7 . 



cbs , h v i(v / B7 (L - x x ) ) - 1 


( 11 ) 


+ <TE»S - =. T*) sinh ( /b 7 (L - Xl)) j 

/Bi ' 

Equation (11) can be inverted to solve for Qi in terms of the measured 
value <T p > - Defining, 

QBAR = Xj Qi k rr R = 2xiHi (12) 


which is the heating rate per unit wire diameter, the following working ex- 
pression results: 

QBAR = A 1 <T p > - A 2 (TEND - TA) (13) 


The constants, Bi, Q ls C 2 , and A 2 are written out in the Appendix. Thus 
QBAR is determined algebraically from the experimental and physical data. 

In the event that Tp(x)/TA is not small compared with unity, the non-linear j 
equation must be solved. A computer program has been written to calculate the 
beam power per unit wire diameter, QBAR, corresponding to a given average 
incremental temperature. This is accomplished by assuming an upper bound on 
QBAR and iteratively converging to the actual QBAR which corresponds to the 
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given average incremental temperature. Comparison of results from the closed 
form solution and the iterative method showed good agreement when a linear 
case was tested. 

2.3 Results 

The incremental change in resistance of the platinum wire was determined 
for various y-positions in the beam (see Figure 4), with set values of target 
region pressure, beam energy and ion current at the beam stop. Via equation (1) 
these values were converted to average incremental temperature rise. Typical 
curves of <Tp> versus arbitrary y-position are presented in Figure 5. The 
beam boundaries are well defined; however, the beam power density distribution 
in the vertical direction through the beam is skewed upward; such was true in 
most but not all beam profiles. This behavior relates closely to beam forming 
and focusing properties determined outside the target region of the system. 

A non-zero "smear" condition is detected at the top edge of the neutral beam, 
which is attributed to straggling ions and neutral atoms. These particles are 
diverted from the principal direction of neutral beam flow at the position 
where vertical deflection plates send the ion beam along the tube above the 
neutral beam to the ion beam stop. Since the neutral beam edges were generally 
distinct and the neutral particle content in the "smear" region was clearly 
small, extrapolation to establish beam diameter, BDIA, was performed on the 
curves on Figure 5. Note that the data points were joined with a smooth curve, 
and a circular beam cross section was assumed, consistent with the mechanism 
of electrostatic ion beam focusing. Besides beam diameter, information 

regarding temperatures TA and TEND wwas'e' needed before employing the data 

r ~yi : :z'z- •' i( 

reduction scheme outlined in section 2.2. The wire end temperature was taken 
to be that of the massive stainless steel mounts shown in Figure 3 to which 
was attached an iron-constantan thermocouple junction. The enclosure tempera- 
ture, TA, was taken as that of the ambient atmosphere in contact with the tar- 
get chamber outer surface. 

Calculations were consequently made to determine QBAR versus y, the 
vertical position of the detector wire in the beam. The curves of Figure 5 
are consequently reduced to those shown in Figure 6. Recalling that QBAR is 
the heat rate per unit wire diameter, it is clear that the total neutral beam 
power is determined by integration after connecting the data points using a . 
jsmooth .curve; integration of the asymmetric beam prof i 1 es was performed Iriechani - 
cally. Neutral beam powers for the cases shown in Figure 6 are. 
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Typical <Tp> data for neutral beams. 





Figure 6. Tower density profiles. 




40.4 milliwatts 

O' - A . 20.3 milliwats 

[^compared with jonbeam powers of 34 ya x 60 kV = 2. 04 watts , 20 ya x 60 kV - 1.20 \ 
watts respectively. The ion beam currents were estimated from measurements 
made with a Faraday cup shielded against secondary electron loss. 


2.3.1 Pressure Effect . Two cases were performed under the following 


conditions: 


Pressure(torr) 

(\>2xib~ e 



Ion current (ya) 

34 

34 

beam potential ‘(kV.)) 

60 

60 

TEND (°R) 

524.7 

530.2 

TA (°R) 

533.7 

534.7 

BDIA (cm) 

1.15 

.74 

f 

neutral beam power (mw) 

13.5 

40.4 


The 9x1 0“ 6 torr case is displayed as the > □ -data in Fi gures 55 earn'd 66 . 'The higher 
pressure data would be expected to yield a greater neutral beam power in direct 
ratio to the pressures ,' i.®., ,4.5 to 1.0. The neutral beam power ratio found is 
in ratio 3.0 to 1.0. This discrepancy is due mainly to the inequality in TEND 
values (a difference of 5.5°R) resulting in a higher neutral beam power for the 
low pressure case. It is not clear why the reference temperature of the wire 
mount, /which did not vary widely while data was taken, should be lower in the 
low pressure case. The BDIA difference is due to focusing irregularities, but 
this should have no effect on neutral beam production so long as the ion cur- 
rents are equal and the beam diameters remain large compared with mean spacing 
between background gas atoms. 


2.3.2 Self-Consistency Test . The wire detector was scanned through both 
the neutral beam and the total beam (neutrals plus ions) in a sequence of ex- 
periments. A separate measurement of the ion beam current was made with the 
shielded Faraday cup. A power balance was attempted such that, 



r 6 + 6. 


(14) 


The symbols represent heating rate due to total beam, neutral beam, and the ion 
beam respectively. The factor g(= 2.63) is a geometric factor gauging that 
length of the beam tube over which fast neutrals were produced by charge 
transfer. For a 40ya, 60 keV ion beam at 2.x 10" 6 torr the results" are, j \ 
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.981 = (2.63) (.018) + 2.400 [watts] 

* 2.447 

The discrepancy is attributed to (1) secondary electron emission problems in 
measuring ion beam current and (2) secondary processes occurring at wire surface 
while scanning the total beam. Ion scattering, secondary emission and sput- 
tering processes on the platinum wire would cause a portion of the total beam 
energy not to be deposited, hence the measured value, Q . = .981 watts, would 
be less than actual total beam power. It is difficult to believe that surface 
processes alone could account for the factor of 2.5 difference in deposited 
power. Both types of errors must share the discrepancy. It is possible that 
there is a spread in ion beam energy, although small for a magnetically analyzed 
beam, such that every ion does not carry 60}' l^.eV?< kd.net.ice energy., ashhascbee"n. as- 
sumed in calculating . In considering the fraction of ion beam which is 
neutralized, the ratio is .018 watts/2.400 watts ~ .8%, a value not inconsistent 
with the V/o reported by McGinty, et al. 3 using a 34pa, 60?vkey..phps.phorust.beam 
with BDIA = .8cm. 

3;m JSPUTTERING MECHANICS 

./• 

A vast literature exists describing ingenious experiments elucidating 
sputtering phenomena induced by primary ions ranging in energy from threshold 
values of a few tens of eVi to the meV„ region (see summaries by Kaminsky, 
and Carter and Colligon). No significant progress was made in this present 
study with respect to quantitative details of argon ion or dopant ion sputtering 
on semiconductors. 

3.1 Semiconductor Sputtering 

A few papers were found (published between 1967 and 1971) mainly addres- 
sing the question of compositional carry-over to substrate and experimental 
sputtering yield comparisons with various theories. A paper by Comas and 
Cooper reports experiments where single crystal GaAs was sputtered by normally 
incident low energy argon ions; composition of sputtered material and yield 
values are offered as results. In section 7.9 of the work by Carter and 
Colligon it is reported that the composition of sputtered material from GaSb and 
InSb samples is basically stoichiometric. In the paper by Sommerfeldt, et al. 3 
30 keV argon ions bombarded a single crystal of germanium. By controlling 
target temperature they were able to study sputter yield and secondary electron 
emission for ordered and amorphous states of the same target. The theory of 





focused collision events along close-packed atomic rows (fbtuis6'n‘' , meq|?a^i ian^f ' • 
did not explain their experimental findings. Another paper taking issue. with 
the focuson mechanism was published by Cowell and Smith. Two other papers, 
relating to semiconductor sputtering, by Hanak also E1;remenkova',>et aZ., are 
listed in the reference section. 

Two conclusions can be drawn from these few recent works: (1) sputter 
theory and experiment have not yet been brought into good agreement (Kaminsky, 
Carter arid Colligon, Sommerfeldt, et al. 3 Cowell and Smith), and (2) composi- 
tion of certain semiconductor materials carries over essentially in stoichio- 
metric ratio (Carter and Colligon, Comas and Cooper). 

3.2 Isotropic Sputtering Model 

If it is assumed that sputtered molecules have no preferred direction of 
departure, as is often the case for polycrystalline materials, the directional 
concentration of sputtered molecules would follow the posing law: 


'T . • 


t 

i ; 




1 

$2 




S * A cos e 

where S n is the number of molecules leaving in the direction e, A is a 
constant and e is the angle between the surface normal and the molecular 
Aside from a difference in velocities, this type of scattering is exactly ' 
analogous to gray diffuse radiation theory. As a result, the sputtered yield 
between two finite areas can be determined using radiation enclosure theory. 

According to Hamilton and Morgan, the fraction of molecules leaving^ 
differential area dA x and arriving at an arbitrarily inclined rectangular, i 
area A 2 (see Figure 7) would be eV. . 


•4*VJ* . V. 

V." 


", 

v 




i 


dh^h 2 ~ 27(dATT 


'tan -1 M +[• — H c° s Lz 1 - V an" 1 


Vl + H Z f-v 2 H COS 
\ M cos <t> 


* + H z - 2H cos <pj + sin z <j> 

+ tan" 1 ( cos — V ) 

\M Z + -sin z <j>y j 


tan' 


H - cos <j> N 
1 Z + si n z <j>. 


where M = m /a and = "/a, (m, n and a are shown in Figure 7). Conversely, 
the number of molecules leaving area A 2 and arriving at dAj would be: 


(16) 


n A? F, 


dAi 
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where n is the number flux leaving A 2 


-r; ^ ^ 

second). 

From radiation theory, 


F '-,4Ai p 

A 2 + dA x Qfep 4&i. * A 2 

- <r" 

Hence, 



(molecules per square cm per 

( 17 ) 


nA 2 F A 


2 


-> dAj 


where 


= ^r 

2tt [ 


tan 


-»M + H: cos 1 tan' 1 r./!L£-<C« 


‘t 1 Sf-i- (tan- 1 + tarr 1 ^±) 


t. 

1 ' 


(18) 


Y 2 = (1 *';H 2 - 2 H cos ij>) and B 2 = M 2 + sin 2 £ 

Consequently, the sputtering yield on a surface can be estimated by calculating 


F»>>* , n for-a-U-- dAr on the receiver -surface. -' 

A^ -*■ dAi ? • • a,A+ ufi T* 0 " — 


If more than one target material is used simultaneously and each material 
is a separate rectangular area, the relative concentrations can be determined 
by breaking the target up into an appropriate number of rectangular areas and 
analyzing each rectangle separately. Then, the. concentration at each dA x 
can be determined for each constituent and the relative concentrations can be 
determined by dividing the number of molecules of a particular constituent 
by the total number of molecules. 


As a test case a target-substrate system as shown in Figure 8a; was 
studied. Initially, the target was assumed homogeneous and the local yield 
was determined by dividing the substrate (receiver) surface into 0.25 x 0.25 cm 
differential areas, then calculating F^ 2 ^ dA to the, center of each "dif- 
ferential 11 area. The relative yields across the substrate surface are shown 
in Figure 8b, s where S/S^gg^ i s the fraction of sputtered material arriving 
at the substrate. 

As a second case, it was assumed that the target was divided into two 
equal areas -- a left half and a right half -- of dimension -j x T w * It 
was assumed that the target materials A and B each had" the same yield? 
under ion bombardment, hence the flux expressions (17) and (18) are valid. 
Naturally, the same total concentration was* computed as before. The relative 
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concentration of A -- the left material -- on the substrate was calculated 
via a computer program. Reasonable results were obtained in tabular form. 

Since no experimental results were available to compare with model calculations, 
the tabular values are not presented here. 

4. r CONCLUSIONS 

The main contribution of this work was to establish a method of determining 

t 

neutral beam flux, dose, or power density. An energy deposition technique was 
used, and chord-wise sampling of beam power density distributions gave some 
fine structure information. Certain secondary processes at the detector wire 
surface limit the accuracy of the technique. The "resolution" of the detector 
was not determined, but neutral beams of approximately 1% ion beam equivalent 
current were detected and beam profiles measured. Integration of beam profiles 
yielded total neutral beam power. 

The method could be further studied in order to establish the accuracy 
of the technique. For the instrument to be used as an on-line monitoring 
device, an automatic scanning and resistance read-out system could be devised. 

As conceived the instrument serves mainly as a research tool. 

A model successfully representing sputter material transport from target 
to substrate was constructed. -Some analytic results are presented, but no 
experimental data was found to compare with model calculations. 
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APPENDIX 


Constants involved in the linear solution of equation (2) for ultimate 
determination of r .QBAR are listed here. 


8 e a (TA) 3 . 

Rk 


Bi 

Qi 

Fi = ,(;(Qi/B 1 ) - (TENDA - TA) [sinh (/B7(L - Xj)) 


Oj 

k-rrR 


cosh (/bTxi ) 
sinh (/Mi) 


+ cosh (/B7(L - x x ))]>/{ sinh (/B7(L - x x )) 


+ cosh (v^b 7(L - xj)) cosh (/^x.jJ/sinh (/Mi)M 


- 4 ' 


C 2 = { -Fi cosh (v^ 7(L - Xl )) 


- (TEND - TA) • sinh (/B7(L - XjJJj/sinh (/Ml) 

A x = (k tt R L)D 2 

A 2 = k tt R D 2 (i^)[sinh (/B7(L - x : ) • -- s - h - ^— x - x ^ 
v sinh (/Ml) 

+ cosh (/b 7 (L - Xx))]/Dx 

Dx = (sinh (/BV(L - xx)) + cosh (^(L - x x )) • cosh 

{ sinh (/Mi) 
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